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The  coefficients  of  the  thermal  conductivity  (* )  and  first  viscosity 
(f?)  in  thin  helium  films  are  evaluated  explicitly  as  a  function  of  tempera¬ 
ture  via  phonon-phonon,  phonon-roton  and  roton-roton  scatterings.  Above 
about  0.8  K,  phonon-roton  scatterings  and  five-phonon  processes  are  the 
main  contributors  to  both  coefficients.  Below  about  0.8  K  both  coefficients 
increase  exponentially  with  decreasing  temperature.  At  temperatures 
below  0.3  K,  has  a  T“s  dependence,  while  shows  exponential  and  T 
dependences.  In  the  case  of  the  former  is  due  to  phonon-roton  scat¬ 

tering  and  the  latter  originates  from  three  phonon  processes.  The  coef¬ 
ficient  k  from  roton-roton  scattering  varies  as  T  1  ,  and  the  roton  part 
of  the  first  viscosity  is  independent  of  temperature. 


1 


PACS  NO.  05.60  +  w,  67.40  Pm,  67.40  -  w. 

*  To  whom  correspondence  should  be  addressed. 


I.  Introduction 


Since  the  works  of  Landau  and  Khalatnikov  on  kinetic  phenomena  in 
superfluid  helium,  there  has  been  continuous  interest  in  thermal  conductivity 
and  viscosity  for  bulk  liquid  4  He  and  3  He  at  low  momenta  and  at  low  tempera- 
tures.  Recently  Kirkpatrick  and  Dorfman"*  obtained  transport  coefficients 
for  very  low  temperatures  (naXJ  »  I)  and  for  moderately  low  temperatures 
(naXJ  «  1)  on  the  basis  of  their  kinetic  equations  for  a  dilute  superfluid, 
where  n,  a  and  X  represent  the  number  density,  the  s-wave  scattering  length 
and  the  thermal  wave  length,  respectively. 

In  the  case  of  thin  helium  films,  the  dissipation  of  superfluid  flow  and 

4) 

thermal  conductance  have  been  investigated  by  many  authors.  Ambegaokar 

5) 

et  al.  have  predicted  that  the  effective  thermal  conductance  in  thin  helium 
films  has  exponential  dependence  on  temperature  for  T  <  and  diverges 
exponentially  for  T  ,  where  is  the  thermodynamic  Kosterlitz-Thouless 

temperature.  More  recently  the  depairing^  and  the  depining  ^  vortices  which 
give  the  power  law  and  the  exponential  dependence  of  thermal  conductance  have 
been  investigated  experimentally  by  Gasparini  et  al. ,  who  confirmed  an 
exponential  dependence  rather  than  the  power  law  for  T  <  T  .  However,  there 
is  much  less  information  about  thermal  conductivity  and  viscosity  in  thin 
helium  films  at  low  momenta  and  at  very  low  and  moderately  low  temperatures. 
What  is  more  important  is  that  a  wrong  normal  dispersion  relation  was  used 
in  Landau  and  Khalatnikov's  well-known  results,^  while  the  correct  disper- 
sion  is  anomalous  dispersion.  For  this  reason  and  in  view  of  recent 
experimental  development  on  helium  films,  we  present  in  this  paper  new 
results  on  the  thermal  conductivity  and  viscosity  of  superfluid  helium  films 
through  the  theory  of  kinetic  phenomena  developed  by  Landau  and  Khalatnikov. 


We  evaluate  thermal  conductivity  and  viscosity  within  three  ranges  or  tem¬ 
peratures:  0.3K  <  T,  0.3  K  <  T  <  0.8  K  and  T  >  0.8  K.  In  these  temperature 
ranges,  scattering  depends  on  the  nature  of  interactions  of  elementary  excita¬ 
tions,  i.e.,  phonon-phonon,  phonon-roton  and  roton-roton  interactions.  The 
scatterings  which  govern  the  transport  processes  and  kinetic  coefficients  of 
thermal  conductivity  and  viscosity  can  be  determined  by  the  characteristic 
time  r  of  scattering. 

In  the  present  paper  we  shall  treat  a  thin  helium  film  as  two-dimensional  (2-D) 

—  less  than  three  atomic  layers,  namely  one  statistical  layer  of  3.6  X  —  and  neglect 
substrate  effects.  In  the  calculations  we  shall  use  the  2-D  excitation  dispersion 
relation  obtained  microscopically : 

e(p)  -  c„p[l  +  5,  pJ  -  <5,  p4  +  ...  ]  ,  (1.1) 

e(p)  -  a  +-ij-  (P-P0)J,  (1.2) 

where  c,,  is  the  sound  velocity,  and  A,  fi  and  P0  are  the  roton  parameters.  To 
evaluate  the  coefficients  of  thermal  conductivity  and  viscosity,  we  shall  first 
calculate  the  scattering  cross  section  for  the  various  interactions  in  Sec. 

I.  Then  we  shall  evaluate  the  characteristic  times  corresponding  to  the 
various  scatterings  by  solving  the  collision  integral,  and  then  obtain  ther¬ 
mal  conductivity  in  Sec.  HI  and  viscosity  in  Sec.  IV.  Finally,  we  shall  give 
results  and  discussion  in  Sec.  V  in  terms  of  graphs  and  tables. 

II.  Scattering  cross  sections  and  differential  decay  rates 

In  this  section  we  consider  the  interactions  of  elementary  excitations 
by  the  second  quantization  method. To  obtain  characteristic  times  corres¬ 
ponding  to  the  three  interactions,  we  first  evaluate  the  scattering  cross 


.'■iVll.tt'Lf. 


section  or  the  differential  decay  rate,  which  are  directly  related  to  the 
collision  integral.  The  collision  process  of  phonon-phonon  interactions 
includes  three-phonon  process  (3PP,  Pj  =*  p2  +  p3 ) ,  four-phonon  process  (4PP, 
Pi  +  P5  ^  P3  +  p4 )  and  five-phonon  process  (3PP,  p,  +  p,  51  p3  +  p4  +  p5  ) . 

The  differential  cross  section  and  differential  decay  rate  for  3PP  in 
two  dimensions  are  defined  as 


d„  .  (  »  )  ,  <FIK,  I  I>|  -  5(Ef-  Et)  T5fyr  , 


<  F  I  H  I  I 


5 (E_  -  E. ) 


(2ffh): 


(2.1) 


(2.2) 


The  3PP  consists  of  two  processes:  the  direct  process  of  emission  of  a  phonon 
P3  by  Pi  =  Pi  +  p3  and  the  reverse  process  of  absorption  of  a  phonon  p3  by  p, 

+  p3  =  p,  .  The  transition  amplitude  between  the  initial  state  I  I  >  and  the 
final  state  I  F  >  is  given  by 


<  F  !  K,  !  1  >  =  5  (P)  -  p  -  P,  )  (  (  5  (P,  *P,  ) 

f  90  3/2  ^0  n  r 3 


+  5  ^  Jp  (  H  Pi  PjPjl^Hn^  (n_^  +  1 )  (n_^  +1)}  5, 

®  ®  n  n 


P|  Pj 


(2.3) 


<  F  '  ^  '  I  >  =  -y  ~ y  5  (P,  -  P,  -  PS H (  )  **  (P,  *Pj > 

4  (2S)  p°  pi  pJ 


+  3  ^  H  ^  p,  p,  P,  )4H(n.  +  On.  ^  , 

Pi  Pi  Pj 


(2.4) 


where  n^  is  the  distribution  function  of  phonons  with  momentum  p.  Then  the 


total  decay  rate  in  both  processes  becomes 


I 


WD  =  (u+1)2  PiPaP  n-*  (n-  +  1)(n-+  1)5(£f  "  EI)  (2rrh) 2  ’ 

°  P,  ?!  P3 


WR  =  2 hp. 


(u+l)J  j  p,p,p(n^+  l)n_  n_  5  (Ep  -  E];)  (2nh)‘ 
1  Pi  Pi  Pi 


(2.5) 

(2.6) 


where  u  is  the  Griineisen  constant  given  by  (P0  / ^  )  (^Co  ^po  )  ■ 
In  4PP  the  transition  matrix  elements^ ^  are  given  by 


<  F  I  X  I  I  >  =  <  F  I  X4  I  I  >  +  2 

J 


<  F  I  Jf,  I  J  >  <  Jjjfj  I  > 


EI  -  EJ 


(2.7) 


where  the  main  contribution  is  due  to  the  second  term,  which  becomes  large 
as  a  result  of  vanishing  denominator.  This  corresponds  to  the  case  where 
5,  is  neglected  and  the  scattering  between  phonons  is  collinear.  Therfore 
we  should  not  discard  the  6l-term  in  Eq.  (1.1).  Taking  account  of  5[ p  , 
which  is  much  smaller  than  unity,  and  small-angle  scattering,  we  can  find  the 
maximum  value  for  <  F  I  XI  I>  ,  in  small-angle  scattering  all  the  phonons  are 
moving  in  the  same  direction,  and  the  conservation  of  momentum  and  energy 
flow  hold.  Thus  the  phonons  moving  in  a  given  direction  attain  equilibrium 
with  each  other  much  faster  than  the  phonons  in  other  directions.  The  differ¬ 
ential  cross  section  (Eq.(2.1))  for  4PP  becomes 


do(p,  p,  ,  P.P,  )  -  [ 


(  HL  1  |  <  f  I  Kl  I  >  I  6  <  £+ei  1  )  (2jrh)v  * 
v  he*, 

(2.8) 


Under  the  condition  p  «  p,  we  obtain 


(u-H);  pi  pi'p' 
8  n  h3  P0JCo  P 


i'p'  S  {e(p)+e(p,  )-e(p')-£(Pi')} 


(i  -  coSe  -  36,  p,3  y 


dp'  . 


(2.9) 


As  for  4PP  above,  Che  direction  of  momenta  of  the  colliding  particles 
is  not  changed  in  the  5PP  case.  The  transition  matrix  elements  in  second- 
order  perturbation  are 


I  i  X  i  I  -TC,  I  n  X  g  1  7C,  I  F  > 
<EI  -  EiHEi  -  V 


(2.10) 


where  some  terms  contain  the  vanishing  denominators  under  the  conditions  of 
=0  and  collinear  scattering.  The  5PP  has  the  maximum  probability  in 
small-angle  scattering  and  leads  to  equilibrium  for  the  phonons  moving  in  a 
given  direction.  Rather  than  calculating  Eq.(2.10)  tediously,  we  make  use 
of  the  kinetic  coefficient  given  by  Landau  and  Khalatnikov.  The  change  of 
rate  per  unit  time  in  the  phonon  numbers  is 


{n,  n,  n3  (0,4-1)  (n, +1)  -  (n,  +1 )  (rij  +l)o,  n,  }  dw 


(2ffh)6 

(2.11) 


where  dw  is  the  differential  decay  rate  defined  by  Eq.(2.2),  which  is  pro¬ 
portional  to  p: ,  N  ^  can  be  expressed  in  terms  of  the  kinetic  coefficient 

T  as 
ph 


n,=  - r  m  ,  , 

ph  ph  ph 


(2.12) 


where  M  ,  is  the  chemical  potential  of  phonons.  The  phonon  distribution 
ph 

function  n  has  small  deviation  from  the  equilibrium  distribution  function  r\,  . 
The  distribution  function  n  can  be  expanded  as  a  function  of  chemical  poten¬ 
tial  as  follows 


n  -  no 


no  (no  +1)  . 

V 


(2.13) 


-5- 


To  simplify  Eq.  (2.11),  we  replace  Uj  ,  n,  ,  a,  and  n  by  Che  equilibrium 
distribution  n,0  ,  n^,  ,  n30  and  n^  ,  and  instead  of  n5  we  substitute  Eq. 
(2.13)  into  Eq.  (2.11).  Then  Eq.  (2.11)  becomes 


ph  '  k-gT 


nto  +  1 )  (n*>  +D  dw 


(2ffh)< 


(2.14) 


Since  n40  and  nJ0  are  much  smaller  than  unity,  they  may  be  neglected  without 
any  disturbance  in  Eq.  (2.14),  and  the  integration  of  dw  is  replaced  by  aver¬ 
age  w.  Then  Eq.  (2.14)  becomes 


ph  -  kgT 


nio  n*  ^  w 


'Pi  Qp2  0-P3 
(2nt\  )6 


(2.15) 


Since  w  and  JC  are  proportional  to  ~  p  and  ~  p3fl  ,  respectively,  the  integral 


of  Eq.  (2.15)  is  proportional  to  p8  .  Averaging  the  momentum  of  phonon  gives 


r  .  =  aT  . 

ph 


(2.16) 


Here  a  is  constant,  which  can  be  determined  experimentally  by  the  attenuation 
coefficient  of  ultrasonic  waves. 

Now  we  return  to  the  scattering  of  phonons  by  rotons.  We  consider  a 
roton  in  the  presence  of  the  phonon  field.  We  can  treat  this  roton  as  a 
particle  in  a  moving  liquid  4  He.  Therefore  there  appears  an  additional  term 
-  P*v.  This  can  be  written  in  symmetric  form 


-  j  (P  •  v  +  v  •  P  ) , 


(2.17) 


where  P  and  v  are  the  momentum  and  velocity  operators.  The  phonon  field 


///  r  .  •  ,• 


WWW 


changes  the  density  of  the  medium,  and  thus  we  may  expand  the  roton  energy 
in  terms  of  the  density  p  '  to  second  order  (  p  '  =  p  -  P0 )  as 


K 


ro 


~  \  A 

op 


ir (p  -  p«)!?  { 


A  + 


t <p  - 


p,  >■» 


|2 


+  . 


♦ 


(2.18) 


where  X  ^  is  given  by  Eq.(1.2).  Since  the  value  of  most  rotons  is  close  to 
P0  ,  we  may  neglect  (P  -  P0 )  and  replace  P  by  P0  .  We  may  also  drop  the  term 
(3A /3p  )p  '  ,  which  is  much  smaller  than  (2.17).  Then  the  interaction  energy 
between  phonons  and  rotons  can  be  expressed  as 


V  =  -  -  ( P •  v  +  v*P) 


32A 
dp 3 


12 


(2.19) 


We  note  that  the  terms  in  the  second  bracket  of  Eq.  (2.19)  have  magnitude 

12) 

on  the  order  of  10-1  to  1  in  3-D  liquid  helium 

When  the  roton  cnanges  momentum  P  to  P',  it  absorbs  a  Dhonon  with  momen¬ 
tum  p  and  emits  a  phonon  with  momentum  p'.  In  these  processes  we  may  consi¬ 
der  two  intermediate  processes,  i.e.,  (I)  P  +  p  -*•  P '  =  P  +  p-p'  and  (II) 

P  -  p'  -*  P'  =  P  -  p'  +  p.  Since  the  roton  momentum  is  much  larger  than  that 
of  the  phonon,  we  may  view  this  interaction  as  similar  to  that  between  heavy 
and  light  particles.  The  momentum  and  the  energy  conservation  law  in  colli¬ 
sion  processes  can  be  written  as 


CP  +  (p  -  P0  )5 


cp' 


+  ar  (> p  +  P 


-  p’i  -  p0)j  . 


(2.20) 


Under  the  conditions  p,  p'  «  P0  and  e  *  cp  «  3  pc2  ,  Eq.  (2.20)  becomes 


-n'  = 


P-P 


2mcP0 


'.P0m-(pn  -  p'n’))  =  {m(n-n')} 


(2.21) 


where  m,  n  and  n'  are  unit  vectors  directed  along  P0 ,  p  and  p',  respectively. 
Therefore  energy  conservation  implies  p  s  p'.  This  means  that  the  light 
particles  do  not  change  the  magnitude  of  momentum  but  change  its  direction. 
Taking  account  of  p,  p'  «  P0  and  P  -  P0 ,  the  matrix  element  (Eq.  (2.1))  in 
second-order  perturbation  can  be  obtained  as 


V'tf  nu 


■n)  +  (m*n')]  (n*n’)  -* — (m#n)‘  (ra*n')2 

Me 


(2.22) 


_£1  r  2lL+  k  {  il2_):i 

P0  c  (  3p2  M  ^  dP  '  ‘ 


(2.23) 


Substituting  Eqs. (2.22)  -  (2.23)  into  Eq.  (2.1)  and  performing  the  inte¬ 
gration  over  p',  we  obtain 


da 


L2_E_ 


8ffh3pJ  c2 


{[(m*n)  +  (m*n' )]  (n'nf )  +—  (m*n)2  (m*n'  )2 

Me 


+  A}Jd0  . 
(2.24) 


Averaging  Eq.  (2.24)  over  all  directions  of  roton  momentum,  we  finally  get 


do  =  teSp/ci  {  (1  +  co^)cob>  +  T55-  (  ( 


128  ^  mc 


35  cos4  0  +  3sin  & 


+  30  cos2//  sin2 1  |  ^  (3  cos2'/'  +  sin2'//)  +  A2  }  d*  ,  (2.25) 


where  is  the  angle  between  the  incident  and  scattered  phonons  (see  Fig.  1). 

We  shall  now  examine  roton-roton  scattering.  Since  the  character  of 

interaction  between  rotons  is  not  known,  we  may  assume  the  short-range 

13) 

potential  given  by  Landau  and  Khalatnikov  and  take  the  roton  interaction 
to  be  a  delta  function  potential 


-8- 


We  construct  the  symmetrized  pairwise  plane  waves  over  incoming  and  outgoing 


rotons  as 


*(P,  ?,  )  = 


L  j  ^ (P-r+P,  t,  )  |(P-r,+P1*r) 


v2  1 


w(p* ,?;> 


i  i  n 

v2  T  e 


i(?  '  -  V  *?.'  *r,  ;  i(P’.r,  +?,'  -r)  i 


The  differential  decay  rate  from  before  to  after  collision  is 


2n  . 

dw  =  —  I  V 

h  AF 


2  d? ' d? ' 

I  «(E  +  E,  -  E*  -  E,’)  , 


where  the  matrix  element  V  is 

AF 


VAF  =  Vo  S  j  (P>Pi  }  5  (r_r<  )^( P',P;)drdri 


_u  f  ^(P'+P.'-P-P,  )-r  _ 

2  V  S  ^  e  dr  . 

0  i 


Performing  integration  over  P'  in  Eq.  (2.28),  we  obtain 


i  071  . 

dw  *  -r—  r  v  I 


1  5  (E  +  E,  -  E*  -  E')  — lE.'— 

1  1  ’  (2?rfi)J  • 


.  (2.27) 


(2.28) 


(2.29) 


(2.30) 


-9- 


i 

); 

i. 


Dividing  Eq.  (2.30)  by  che  relative  velocity  of  rotons  given  as 


3_E  _  3Ej_  |  (2.31) 

—  —  » 
dp  d  p, 


we  can  obtain  the  differential  scattering  cross  section. 

As  mentioned  earlier  in  the  discussion  of  the  scattering  of  phonons 

by  rotons,  most  rotons  have  the  value  close  to  P0  .  Thus  the  change  of  momen¬ 
tum  after  collision  is  very  small  in  comparison  with  P0 .  Let  us  take  8  as 
an  angle  between  the  incident  rotons  with  momenta  P  and  P,  before  collision 
and  introduce  the  variable  f  (Fig.  2).  Then  the  roton  momenta  after  collision 
can  be  expressed  as 


e  e 

P'  =  P0  +  fx  cos  -  +  fy  sin  -  , 

B  B 

P.  =  P0  -  f  cos  -  +  f  sin  ■=. 

1  0  x  2  y  2 


(2.32) 


„  ,  ,  ,nri  I  "fl  «  P  .  From  the  conservation 

Here,  we  have  made  use  of  P  -  F;  -  r0  ana  i  ri  ^  r0  . 

of  energy  we  have 


f  2 
x 


2  B 

cos  2 


+  f  2 

y 


Sin  r 


i  (P  -  Po):  +  5  <pi  -  p0)3 


(2.33) 


To  obtain  the  total  scattering  cross  section,  we  integrate  Eq.  (2.30) 
together  with  Eqs.  (2.31)  and  (2.32)  over  the  momentum  of  the  scattered 
particles  to  get  the  total  scattering  cross  section 


a 


Afi  I  V  I  2 
or 


— I h3  sin# 
9  P, 


(2.34) 


(2.35) 


and  the  average  collision  cime  between  rotons  becomes 


1  Jr  sr  4/iiV  I2 

-L  =  ff|il_i^iN=  - or_ 

Cr  bp  bp,  '  r  h3sinfl  ‘r  * 


where  N  is  number  of  rotons  per  unit  area  given  by 


B  I  P0 


-a/k  T 


[  2jt  j  h2 


(2.36) 


HI.  Coefficient  of  thermal  conductivity 


The  equilibrium  distribution  function  r\j  of  excitation  satisfies  th’ 


kinetic  equation 


bn  ,  bn  b3f  bn  b-K  T,  . 

bt  br  bp  bp  br 


(3.1) 


with  vanishing  collision  integral.  When  the  equilibrium  is  disturbed,  we 


assume  that  the  nonequilibrium  distribution  function  n  deviates  slightly 


from  equilibrium.  The  small  deviation  can  be  determined  by  the  first 


derivatives  of  the  velocities  v  ,  v  and  the  thermodynamic  variables,  since 

n  s 


function  can  be  written  as  n  =  n,  +  n,  ,  where  n,  «  .  Putting  n  into  the 


left-hand  side  in  Eq.  (3.1),  it  is  sufficient  to  keep  only  the  differentiation 


>f  Oq  since  the  derivative  of  n,  makes  higher  derivatives  which  can  be 


neglected.  For  the  collision  integral  on  the  right-hand  side  we  only  keep 


the  linear  terms  in  n,  .  With  the  help  of  the  continuity  equation,  entropy 


'■  -»  JV  "»  V.V  ’ 


IT*  V  * 


-  rys 


T*  ^  WJBWA  77^77  VHWV 


equation  and  superfluid  equation  of  motion,  we  can  write  the  kinetic  equation 
(Eq.  (3.1))  as 


n'  r  ,  1  5T  .  9e  -  •* 

T  9p"  £  ”  9p"  ^  V*<J  ~ 


<j  -*v  +[^Bp  +Hs)£-lf  i^-vn 


+  Y‘*(p^-e^:)-^§  V(p.vn)  }  =  J(n,  )  , 
n  3  p  9  p 


9  £ 


9  e 


n 


(3.2) 


where  n'  =  -  n(n  +  1). 

When  there  exists  a  temperature  gradient  in  superfluid  4  He,  there  is  not 
only  transport  of  heat  but  also  an  irreversible  heat  flow  which  can  be 
expressed  by  the  coefficient  of  thermal  conduction  given  as 


q  =  -  k  7T . 


(3.3) 


Comparison  of  Eq.  (3.3)  with  Eq.  (3.2)  for  thermal  conductivity  k  gives  the 
kinetic  equation  with  a  temperature  gradient, 


— Mpfi- 


ST 

P 


n 


9  £ 
9  P 


=  J(nj  ), 


(3.4) 


L 


where  9  is  the  angle  between  p  and  VT.  The  phenomena  associated  with  ther¬ 
mal  conduction  in  bulk  liquid  4  He  have  aspects  in  common  with  thermal  trans¬ 
port  properties  of  ordinary  classical  liquids.  However,  there  are  specific 
features  of  thermal  transport  which  are  connected  with  the  unusual  elementary 
excitations  of  liquid  4  He.  The  left-hand  side  of  Eq.  (3.4)  is  always  zero 
for  a  pure  phonon  gas,  and  thus  the  corresponding  thermal  conductivity  vani¬ 
shes. 

The  thermal  conductivity  *  consists  of  two  parts,  i.e.,  the  phonon  k 


ph 


and  the  roton  k  : 

r 
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Let  us  first  consider  the  roton  part.  This  is  determined  by  the  roton-roton 
scattering  process.  As  mentioned  earlier  in  Sec.  2,  the  character  of  the 
roton-roton  interaction  is  not  well  known,  and  we  thus  obtained  the  average 
collision  time  t  under  the  assumption  of  the  delta  function  potential. 

Since  we  have  only  to  know  the  temperature  dependence  on  we  replace  the 
collision  integral  by  the  following: 


J  (n) 


(3.6) 


Substitution  of  Eq.  (3.6)  in  Eq.  (3.4)  yields 


n  ~  =  - 


E__yT.[p|I_£i£]  tr. 


V 


(3.7) 


ap 


Substituting  Eq.  (3.7)  into  the  expression  for  the  energy  flow 


£(P)(n  -n,)  dP 
3P 


(2rrh)2  ’ 


(3.8) 


and  comparing  this  result  with  Eq.  (3.3),  we  get 


1  f  .  ST  3e  N 

k  =•  t_  —  I  n*  e  ) 


r  '  r  ZkgT2 


dP 


ap 


p  -  '  (2nh  )J 

n  °p 


(3.9) 


With  the  help  of  n1  =  -  n  and  Eq.  (1.2),  we  finally  obtain  the  roton  part  of 
thermal  conductivity  as 


«r(T) 


r  u 


2mT 


1  + 


o  f  a_  'i  r  a  ,  ^  r 

T  t~)  -p~F  +  2 

n  B 


J  'j  ' 

(3.10) 


According  to  the  scattering  processes  in  Sec.  H,  the  collision  integral 
J (n)  becomes 


J(n)  =  J3pp(n)  +  J4pp(n)  +  J5pp(n)  +  J^n) .  (3.11) 

The  4PP  do  not  change  the  total  number  of  phonons  but  have  a  characteristic 
temperature  T'  in  a  given  direction,  which  is  different  from  the  temperature 
T  in  the  equilibrium  state.  The  law  of  energy  conservation  yields 

p 

J4pp(n)  pdp  =  0  .  (3.12) 

The  total  number  of  phonons  traveling  in  a  given  direction  is  changed  by 
small-angle  3PP  and  5PP.  Therefore  the  distribution  function,  which  not  only 
depends  on  temperature  T'  but  also  the  chemical  potential  a',  can  be  written 
as 


n  =  exp  [(<*'  +  pc/k^T’)  -  1]  ”  . 


(3.13) 


Expanding  Eq.  (3.13)  as  a  function  of  T'  -  T,  we  can  express  Eq.  (3.13)  in 
terms  of  equilibrium  distribution  function 


5n  =  n  -  no  =  -no(no  +1)  {a' 


PC  T'-T  | 

V  T 


(3.14) 


Since  the  left-hand  side  of  Eq.  (3.4)  is  involved  in  9  ,  a'  and  T’-  T 


depend  naturally  on  d  .  To  solve  Eq.  (3.4)  we  should  take  the  forms 


a  1  -  a  cos  d  , 


(T'  -  I)/T  =  P  cos  6  , 


(3.15) 


where  a  and  (3  are  constants  to  be  determined  by  the  kinetic  equation.  Con¬ 
sidering  the  conservation  of  phonon  numbers  in  a  given  direction  and  conserva¬ 
tion  of  energy,  the  integrations 


J3p^n)epdp, 


vanish  and  Eq.(3.4)  becomes 


J4pp(n)£pdp, 


J5PP(n)£pdp 


n' 

3T  r 

ST 

-  3e  i 

J  kBT 

COS  T  [ 

fix 

P-' 

n 

"  3p  ‘ 

f  n' 

cos  0  — 
ox 

r  ST 

d£ 

1  V 

fpp- 

n 

"  L  3P 

)  £pdp 


(n)£pdp. 


j  k^T - 3x  pn  3P  J  Ph 

From  Eqs.(2.5)  -  (2.6)  the  collision  integral  for  the  3PP  becomes 


(3.16) 

(3.17) 


31'jp  ( n ) 


3*7  (u  +  D2  ,  p,p5p  6n(nPi;>  -  npJ6  (£p  -  yyrffr 

(3.18) 


where  6 n  is  equal  to  (n  -  n  )  and  n  represents  the  equilibrium  distribu- 

P  Po  Po 

tion  functions  for  the  phonons  with  momentum  p.  Making  use  of  Eqs.  (2.12)  - 
(2.13),  the  collision  integral  for  the  5PP  can  be  written  as 


J5pp(n)pdp  =  In  h1  kQT  a T 


(3.19) 


We  now  evaluate  the  collision  integral  in  the  scattering  of  phonons  by 
rotons.  When  a  phonon  with  momentum  p  changes  to  momentum  p'  directed  at 
angle  6  '  after  collision  by  a  roton,  the  probability  per  unit  length  that  a 
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/  -  /..wmv.v.vw 


Hi 


m 


particle  undergoes  collision  is  N  da,  where  N  is  the  roton  distribution 

r  r 

given  by  Eq.  (2.36)  and  da  is  given  by  Eq.  (2.25).  Then  the  collision 

integral  J  (n)  can  be  written  as 
ph-r 


Jph^(tl)  =  '  NrC°  J  [n(p*  9’  T)  “  n(p’>  9'>  T,)  1  d°' 


(3.20) 


Here,  9  is  the  angle  of  the  incident  phonon  with  momentum  p  with  respect  to 
the  x-axis.  Using  Eqs.  (_3.14)  -  (3.15)  for  the  difference  between  distribu 
tions,  we  obtain 


V<-(“)  ‘  cose  +  \  aiSpfe- 


x  [  y  +  rrr  ( - ^  A  +  2 A2  ]  . 
4  32  (i  c  n  c 


(3.21) 


Substituting  Eqs.  (3.18)  and  (3.21)  together  with  Eq.  (3.19)  and  n'  =  -  ( 

+  1)  in  Eqs.  (3.16)  and  (3.17),  and  performing  the  momentum  integral,  we 
obtain 


1  il  _ST_  .  2 !r  (2)6!£  (6)  _  _J_  f  _  5!?  (5)  ,  +  _£L 

T  3x  (C  "  P  C  }  3!f  (3)4!f  (4)  "  r_L._l  P4!f(4)J  r  ^ 


Here  r  .  f  b  and  t  are  the  characteristic  times  (see  the  Appendix), 
3PP  5PP  ph-r 


'  U'JU'Vfe'U 


which  characterize  each  collision  processes,  given  by 


Tph-r(T) 


6!f (6)  P°:  (kBT)J  ,  1 
3 !f  (3)  IhW  r  l4 


6!f  (6)  2nh2  c2  a 


■77  (^-)J  +  (  )A  +  2A2]  , 

32  mc  pc  * 


r5pp(T)  3  (3)4  !f  (3)  ^ 


(3.24) 


(3.25) 


1 _  _  2!f (2)6!f  (6)  (u+I)  kE 


r 3pp(T) 


4!f(4)  8ffh3pc 


(3.26) 


Solving  Eqs.  (3.22)  and  (3.23)  for  a  and  P,  we  obtain 


_ 8.071 _ \  iH  JST_ 

0.186  ,  1  ,  1  T  3x  C  "  p  C 

- -  +  - -  +  - -  n 

ph-t  5PP  3PP 


-  A  I  *1  (C  -  -SI.) 

T  3x  ^  P  C}  ’ 
n 


(3.27) 


P  =  -r  + _ U1 1 _  I9T(c__S^) 

ph-r  0.186  ,  1  1  T  3x  V  P  C  1 

-  + -  +  -  n 

^  7  ph-r  T  5PP  T  3PP  - 


-  p  1  ?T  .  ST 

B  T  dx  (C  p  C  )f 
n 


(3.28) 


where  we  introduced  A  and  B  for  simplification.  Substitution  of  Eqs.  (3.27) 
and  (3.28)  in  Eq.  (3.14)  yields 


n  -  «o  =  -  "o  ("0  +  Dcos*  (3.29) 


and  the  energy  flux  (Eq.  (3.8))  for  the  phonon-roton  process  can  be  evaluated 


towuwwiwmiiquuw?^^  v  *>-  *.-  *- iw 


f  3e  s  ap 

j  ^  £  (  n  -  n,  )  cos* 


-4t—  f  C  -  —  )  [  3.290A 
4»rhJC  1  PC'1 


3T 

7.212B  ]  ^ 


(3.30) 


Comparison  of  Eq.  (3.30)  with  Eq.  (3.3)  gives  the  coefficient  of  thermal 


conduction: 


1.803k3  T2 


K  .  (I)  = 


_ 2.310 

T  ph-r  0.186  +  l 

T  ph-r  T5PP 


(3.31) 


We  note  that  the  coefficient  of  the  thermal  conduction  in  bulk  liquid  helium 
is  given  by 


k(T)  =  2  x  io»  i  +  t~3/2  e  '  B  (l  - 

a 


f  „  1  *  °-75 
!  1  +  8  W'5fP 


T  >  0.9  K 


T  <  0.9  K 


(3.32) 


S  N*  V 


IV. 


Viscosity 


In  this  section  we  investigate  the  first  viscosity  through  similar 
calculations  to  those  used  above  for  the  thermal  conductivity.  Equation 
(3.4)  vanishes  for  the  case  of  a  pure  phonon  gas,  and  thus  the  corresponding 
coefficient  k  (T)  becomes  zero.  However,  the  kinetic  equation 


ot 


(4.1) 


does  not  vanish  for  the  pure  phonon  gas  near  zero  temperature.  Therefore, 

we  should  consider  the  contribution  of  the  pure  phonon  gas  to  viscosity. 

Let  us  consider  the  macroscopic  motion  of  liquid  helium  (which  does  not 

— ►  — ► 

depend  explicitly  on  time)  with  velocity  U  and  the  gradient  of  U  directed 
along  the  x-axis.  Then  Eq.  (4.1)  becomes 

v»Vn  =  vx  =  J(n)  •  (4.2) 


The  equilibrium  distribution  function  of  roton  in  liquid  helium  with  velocity 
U  can  be  expressed  bv 


"o 


A  (P-Pp  )2 
kBT  “  2MkB 


(4.3) 


where  we  used  Eq.  (1.2).  Substituting  n  —  and  the  roton  velocity  v  =  (P 
-  P0)/M  in  Eq.  (4.2),  we  get 


no 


Pq  (P-Po  ) 


d  U 

t —  cos  0  sins 
ox 


J3  (n) , 


(4.4) 


where  9  is  the  angle  between  the  roton  momentum  P  and  the  x-axis.  The  colli- 


sion  integral  can  be  replaced  by  Eq.  (3.6)  to  give 


=  -t  (P~P?^P-^  cosd  sind 

r  m  k_T  ^  3x 


n  -  Oq 


(4.5) 


Substituting  Eq.  (4.5)  into  the  following  stress  tensor  of  2-D  liquid  helium 


In  ,  x  dP 

a  =  P  v  (n  -  n. ) 
xy  J  y  x 


(27Th)2  ’ 


(4.6) 


and  performing  the  integral  over  the  momentum  space,  we  arrive  at  the  expres¬ 


sion 


„  _  JY  311 

xy  Cr  8m  ''r  3x  ' 


(4.7) 


Comparing  Eq.  (4.7)  with  the  general  expression  of  viscosity 


au 

°xy  =  \  3x  ’ 


(4.8) 


we  obtain  the  coefficient  of  viscosity  for  the  roton  part  as 


P  2 

n  =  t  n  . 
r  8m  r  r 


(4.9) 


The  equilibrium  distribution  function  of  the  phonon  gas  is 

«.  -!  exp  (  CP  -  H  )  -  1  ]-'  . 

V 

Let  us  assume  the  distribution  function  to  deviate  slightly  from  ,  i.e., 
n  *  +  5n,  which  satisfies  Eq.  (4.1).  Substituting  Eq.  (4.10)  into  Eq. 


(4.2)  we  get 


"o  (no 


4—  cos  9  sin0 
ox 


=  J(n). 


(4.11; 


Solving  Eq.  (4.11)  we  can  obtain  the  phonon  part  of  the  viscosity.  It  is 
necessary  to  consider  J (n)  in  the  various  collision  processes.  Then  the 
collision  integral  J(n)  in  Eq.  (4.11)  can  be  replaced  by  Eq.  (3.10).  Through 
the  same  processes  that  we  argued  in  Sec.  HI,  we  obtain  the  following  two 
equations : 


no  (no+1) 


cos 9  sin0  pdp 
ox 


[J3PP(n) 


J5f>P(n) 


J’  (») 
ph-r 


1  pdp, 


(4.12 


f 


tio  (no  +1) 


cp  5  U 

V  3x 


cos9  sin#  £(p)pdp 


J'  (n)  £ (p)pdp • 
ph-r 


(4.13 


We  can  express  the  distribution  function  by  Eq.  (3.14),  except  that  now  the 
dependence  of  <*'  and  (T'  -  T)/T  on  the  angle  9  is  given  by 


a’  *  a  cos#  sin#  ,  (T ’  -  T)/T  -  P  cos 9  sin#  . 


(4.14 


instead  of  Eq.  (3.15). 

In  Eq.  (4.12)  the  collision  integrals  J^pp(n)  and  J^pp(n)  are  given  by 
Eqs.  (3.18)  and  (3.19)  respectively.  Using  Eqs.  (3.15)  -  (3.20),  we  obtain 
the  collision  integral  J'  (n) : 


,v.vt,v.m  v'.v.r.r.r.rr 


Jph-rU)  =  C0S"  sint*  Nr  Con°  ln<> +  L)  {a  "  (i  'T2f  ) 


—ip 2 f  i  4.  JL  (  Za_  \2  .  5  .  ,  Po_  \ 

8fiJp:CJ  2  32  PC  4  A  PC  ^ 


(4.15) 


Here,  we  note  that  Eq.  (4.15)  is  slightly  different  from  Eq.  (3.21),  because 

Eq.  (4.14)  contains  an  extra  sin#  term.  With  the  use  of  similar  calculations 

for  Eqs.  (3.16)  -  (3.26),  we  obtain  a  0  and  t'"1  as 

ph 


8.071 


0.186  1  1 

r  '  r  t 

ph-r  5PP  3PP 


_ 

1 _  .  dx 


-‘•IS  * 


(4.16) 


!  -  r  '  +  - 

1  phr  0.186 


1.372 


□ » in 

dx  ’ 


(4.17) 


r-  =  M.(6)  P°'kBT2  r  i  ,  _5_  (  Po.)2  +  2  f  h.)  +  2AJ  1 

ph-r  3 !f  (3)  8h3pJc“  1  2  +  32  ^  pC  ’  +  4  ^  pC  ^  2A  1  ’ 

(4.18) 

where  Eq.  (4.18)  is  not  equal  to  Eq.  (3.24).  Substitution  of  Eqs.  (4.16)  -• 
(4.  17)  in  Eq.  (3.14)  yields 


Sn  =  ~  "o  <no  +  D  (A'  -  B'  )  |^  . 

V  9x 


(4.19) 


Combining  Eq.  (4.19)  with  Eq.  (4,16)  and  integrating  over  momentum  space 
gives 


C  r  B 


kDT  , 


xy  16JrhJ  1  C 


—  }  [  3.290A'  -  7.212  B'  ]  • 

C  dx 


(4.20) 


Comparing  Eq.  (4.20)  with  Eq.  (4.8),  we  get  the  phonon  part  of  the  coefficient 


"  -  «  -  j«  1  -  .  .V  .»  V  _r  'ji  •_*  *  A  _m  -  •  , 


fyv.v  s. 


of  che  first  viscosity: 


V  AT)  = 
ph 


7 . 212(kgT; 


2.310 


16ffh2CJ  L  ph-r  0. 186  +  1  +  j _ 

T ph-r  5PP  3PP 


(4.21) 


From  the  theory  and  experiments  of  attenuation  of  ultrasonic  sound  waves, 
we  can  confirm  chat  at  near-zero  temperature  the  contribution  of  3PP  to  the 
viscosity  plays  a  main  role,  and  4PP  establish  only  the  equilibrium  of  energy. 
Therefore,  we  consider  the  3PP  contribution  separately  from  other  contribu¬ 
tions.  The  kinetic  equation  for  3PP  in  Eq.  (4.11)  becomes 


no  (no  +  D  ^  gcosSsinS  =  J3pp(«) 


(4.22) 


Since  the  distribution  function  of  phonons  depends  on  rhe  chemical  poten¬ 
tial  a,  Eqs .  (3.14)  -  (3.15)  become 


5  n  =  n  -  a,,  =  -no(no  +  1)  a ' 


a'  =  a  cos  d  sind . 


(4.23) 


With  a  similar  calculation  for  Eqs.  (3.16)  -  (3.26),  we  obtain  the  stress 


tensor 


[  2!f(2)  ]  2  6!f  (6)  V 


3!4!f (3)f (4) 


16ffh2CJ  3PP  dx 


(4.24) 


From  comparison  of  Eq.  (4.24)  with  Eq.  (4.8)  we  obtain  the  phonon  part  of  the 


coefficient  v  ,  near  zero  temperature  as 

ph 
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1 


(T) 


2.645  ,  3 

7TJ  hJ  C2  ii  '  3PP  * 


(4.25) 


We  note  that  the  coefficient  of  the  first  viscosity  in  bulk  liquid  helium  is 
given  by 


V.  Results  and  discussion 

In  the  previous  sections  we  have  evaluated  tue  scattering  cross  sections 
and  characteristic  times  for  various  interactions  of  the  elementary  excitations. 
Using  these  results  we  have  obtained  the  coefficients  of  thermal  conduction 
and  first  viscosity.  To  investigate  the  temperature  variation  of  the  coeffici¬ 
ents  *  (T)  and  V  (T) ,  we  adopt  the  parameters  which  are  determined  from  the 
specific  heat  data  of  Bretz  et  al.^"^  The  parameters  are  listed  in  Table  I. 
With  this  choice  we  obtained  c  =  84.06  m/s,  which  is  smaller  than  the  value 
157  m/s  of  Hipolito  and  Lobo,^^  but  is  very  close  to  the  experimental  value 

of  (76  ±  2)  m/s  of  Wushburn  et  al.^^  The  parameters  u,  A  and  a  in  Eqs.  (2.9), 

18)  9) 

(2.23)  and  (2.16)  for  two  dimensions  are  assumed  to  be  1.8,  0.425  and 

10  \ 

1.0  x  lO4^  used  by  previous  workers  for  the  bulk  case. 

The  coefficient  of  the  thermal  conduction  is  given  by  Eqs.  (3.10)  and 


-24- 


(3.31).  We  can  confirm  easily  that  Eq.  (3.4)  vanishes  by  considering  only  a 
pure  phonon  gas  because  of  che  peculiar  excitation  of  liquid  helium.  There¬ 
fore,  the  thermal  conduction  depends  on  the  interactions  between  excitations, 
and  it  is  necessary  to  take  into  account  the  phonon-roton  and  roton-roton 
interactions . 

The  roton  part  k  (T)  of  the  thermal  conductivity  is  proportional  to  the 
P 

average  collision  time  t  .  When  9  =  0,  t^  is  zero,  and  for  small-angle  scat¬ 
tering  t  becomes  very  small.  Since  Eq.  (5.10)  is  involved  in  the  unknown 

interaction  potential  constant  V  and  we  have  onlv  to  know  the  magnitude  of 

or 

the  temperature  dependence  for  «  ,  there  is  no  problem  as  to  whether  we  take 

IT 

maximum  t  or  not.  Taking  9  =  jt  / 2  and  the  numerical  parameters  in  Table  I, 
*r(T)  can  be  expressed  as 


n  i q i  x  i n-w  i  i 

*r(T)  =  ~fv~  F"  Y  +  0-221T  +  0.723  +  4.743  (±  +  0.362) 
or  1 


_  4 

1  +  28 . 02 (^  +  )  T'3/2e  T 

_  *.u 

1  +  4632.89  T’7/1  e  T 


(5.1) 


Figure  3  illustrates  «  as  a  function  of  temperature.  The  coefficient  in¬ 
creases  slowly  as  temperature  decreases  and  is  proportional  to  T  1 . 

The  characteristic  times  which  determine  the  coefficient  #  are  8^ven 

by  Eqs .  (3.24)  -  (3.26).  r'*  has  the  temperature  dependence  of  T7/1e  A^kBT> 

ph-r 

which  is  one  power  lower  than  T9n  e  in  bulk  liquid  helium.  This  is  due  to 

sionality,  where  t"1^  has  t*ie  temperature  variation  of  T6  ,  which  is  a  lower 
power  than  that  of  three  dimensions.  The  inverse  of  T^pp  proportional  to 
T4,  which  originates  from  the  anomalous  excitation  spectrum.  Figure  4  is  the 
temperature  variation  of  characteristic  times.  T5pp  comParable  with 


dimen- 


near  about  0.8  K,  and  thus  we  can  conclude  that  above  0.8  K  the  5PP  and 


Tph 

the  scattering  of  phonons  by  rotons  will  mainly  contribute  to  the  thermal 
conduction.  Since  5PP  appear  in  the  inelastic  collision  process  in  the  tem¬ 
perature  range  of  roton-roton  collisions,  we  should  not  take  into  account  this 

contribution  to  «  (T)  below  0.8  K,  and  should  consider  only  the  phonon-roton 

ph 

collision  process.  Therefore,  in  accordance  with  the  regions  of  temperature, 
we  can  express  Eq.  (3.31)  as 


1.803kilJ  CT 

VT>  -  -,w~  <c‘  ' 

^  a 


ph*r 


r  13.419  4-  5.376  1 

1  +  5'376  %  '■ 


T  >  0. 8K 


(5.2) 


and 


V(T)  = 


1 .803koT 

D 

rrh2  O2 


(c  -f) 

n 


ph-r 


=  2.36  x  10'u(l 


_ST  i  t-j/j  g-i"  T  <  0.8K. 

p  c2 

n 


(5.3) 


Since  the  temperature-dependent  term  within  the  parenthesis  of  Eq.  (5.3)  is 
much  smaller  than  unity  for  temperatures  below  about  0.3  K,  it  can  be  neglected. 
Then  we  have 


k  ,  (T)  -  1.093  x  10'7  T-s  , 

ph 


T  <  0.3K. 


(5.4) 


In  Fig.  5  the  coefficient  k  (T)  is  plotted  as  a  function  of  temperature 
based  on  Eqs.  (5.2)  -  (5.4).  As  temperature  decreases  (  <  0.8  K) ,  Eqs.  (5.3) 
plays  a  dominant  role  and  gives  the  exponential  increase.  At  yet  lower  tempera- 


w 


cures  (below  about  0.3  K) ,  Eq.  (5.4)  cakes  pare  in  the  thermal  conducti¬ 
vity,  and  K  ,  increases  rapidlv  and  diverges  as  temperature  tends  to  zero, 
pn 

Comparison  of  Eq.  (5.3)  with  Eq.  (3.32)  reveals  that  the  second  term  in 
Eq.  (3.32)  has  the  same  factor  T-3«  e^/kBT#  This  is  due  to  the  fact  that 
r_lh  and  the  thermodynamic  functions  are  one  power  lower  in  temperature  than 
those  of  three  dimensions,  and  thus  their  ratios  have  the  same  temperature- 
dependent  factor. 

For  viscosity  we  may  give  similar  arguements  as  those  for  thermal  con¬ 
ductivity.  However,  the  main  difference  from  thermal  conduction  is  that  the 
kinetic  equation  (Eq.  (4.1))  does  not  vanish  for  a  pure  phonon  gas  near  zero 
temperature.  Therefore,  we  have  treated  separately  the  contribution  for  the 
pure  phonon  gas  to  viscosity  (Eq.  (4.25)). 

Substitution  of  Eq.  (2.35)  in  Eq.  (4.9)  and  the  choice  of  9  -  rr/2  gives 
the  roton  part  of  the  first  viscosity  as 


h3p  3 

VT)  =  32^r1r 


(5.5) 


which  is  independent  of  temperature.  Since  Eq.  (4.3)  decreases  exponentially 
with  temperature  and  the  mean  free  path  of  roton  increases  at  about  same  rate, 
the  viscosity  becomes  independent  of  temperature. 

At  temperatures  about  above  0.8  K,  the  main  contribution  is  due  to 
5PP  and  the  phonon-roton  collision  process,  and  thus  Eq.  (4.21)  becomes 


/t'i  7.212  .3  , 

%h  T  16*hjq  (kBT)  ph-r 


f  13.419  +  5.376  T’h_r'/T5Pp 


ph-r  [  1  +  5.376  r ph-/T 5PP 


T  >  0.8K  . 


(5.6) 


Wm  PC*  A.  v  ■  *C*  "V  »•  *  ^  .  t  * 


VvOVT v'  OTOv 


At  temperatures  below  0.8  K,  we  consider  only  the  pnonon-roton  collision 
process  and  thus  have 

Vm  ■  (V>‘  T^r 


2.17  x  L0'*  T”1'1 


(5.7) 


For  temperatures  below  about  0.3  K,  we  have  another  temperature  dependence 
given  by  Eq.  (4.25): 


h  (I)  =  5.92  x  10"14  T*1  . 

ph 


(5.8) 


Figure  6  illustrates  the  coefficient  hp^(T)  of  the  first  viscosity  as  a 

function  of  temperature  based  on  Eq.  (5.8).  At  temperatures  below  about  0.8 

K,  i/pk  increases  exponentially  as  temperature  decreases.  However,  as 

temperature  approaches  absolute  zero,  the  roton  density  becomes  zero,  and  thus 

the  contribution  from  the  roton-phonon  scattering  to  viscosity  vanishes  so 

that  the  main  contribution  comes  from  the  3PP.  Therefore,  the  temperature 

U  A/k_T  i  . 

dependence  is  changed  from  T  e  to  I  The  coefficient  increases 

slowly,  and  near  zero  temperature  it  increases  very  rapidly  and  finally 
diverges . 

In  conclusion,  we  remark  that  the  behavior  of  the  coefficients  of  the 
thermal  conductivity  and  first  viscosity  are  very  much  like  that  of  the  bulk 
case.  At  below  0.3  K,  the  contribution  to  the  coefficient  of  first  viscosity 
is  due  to  the  3PP,  which  is  shown  to  have  a  T  ^  dependence. 
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Appendix 


Since  the  calculations  for  Eqs.  (3.24)  —  (3.26)  are  very  similar,  we  shall 
derive  only  Eq.  (3.24).  The  collision  integral  (Eq.  (3.20))  can  be  expressed 


J  (n)  =  N  C  [  (n1  -  no  )  -  (n  -  no  )  ]  da  , 
ph-r  r 


Making  use  of  Eq.  (3.14),  Eq.  (A.l)  becomes 


Jph-r(n)  =  MrC  |  (n«  +  (cos0'  -  cosfl) 


In  Fig.  1  the  relation  between  angles  is  given  by 


cos  O'  -  cosfl  =  cos  0  (  cos  ^  -  1  +  tan#  sin^  ] 


(A.l) 


(A. 2) 


(A. 3) 


Substituting  Eqs. (2.25)and  (A. 3)  in  (A. 2)  and  performing  the  integration  over 


4>,  we  get 


Jph-/n^  =  NrC  COS0  ("0  'jpf  )  °o  ("o  +  1) 


8h3 PJ  C3 


;  +  +ttt  +2a2]  ■ 


The  collision  integral  may  be  written  as 


(A. 4) 


J  (n)  = 

pn-r 


TPh  <P>  ‘ 


(A. 5) 


§ 


b 


ms 


i 


To  obtain  the 


temperature  dependence  of  r  <T), 

ph-r 


integral : 


we  calculate  the  following 


Jph-r(n) 


cpdp/ 


J 


^  rip 


pdp 


(A.  6) 


Performing  the  integrations  over  momentum  space,  we  have 


j  Jph-r(n)  £pdp  =  cos*  \  ^ 

B 


Po 

wy 


2 


6 !f  (6), 


(A. 7) 


:  0 
;  oe 


kRT  3 

Pdp  =  -  3!f  ( 3 )  C  i-S-}  , 


(A. 8) 


where  G  represents  the  bracket  of  the  right-hand  side  in  Eq.  (A. 4),  and  then 
Eq.  (A. 6)  becomes 


cos  6 


N  P 

r 


6!  (6)  p°2y , 

3!  (3)  8h3p2  C2  G’ 


(A.  9) 


Substitution  of  Eqs.  (A. 5)  and  Eq .  (3.14)  in  Eq . 
momentum  space  gives  Eq.  (A.b)  as 

£  =  0  cos0  / r  (T), 

ph-r 


(A. 6)  and  integration  over 


(A. 10) 


Comparing  Eqs.  (A. 10)  and  (A. 9),  we  obtain  the  characteristic  time  r  (T) 

ph-r 


r  (T)  =  \P°  j  1  9  f  Pa.  )  ‘  +  Po  A  2 

Ph-r'  '  3  !f  (3)  8h V  C4  |  4  +  32  <•  7c.  *  +  7c-  +  2A 


(A. 11) 


In  the  case  of  the  first  phonon  viscosity,  we  should  replace  Eq.  (A. 4)  by 


Eq .  (4,15).  Through  the  same  calculation  we  get 


r£r(T)  - 


6!?  (6) 
3 !?  (3) 


r  0  3 

8h3p2  C4 


1  +  A  (SSL 

2  32  '•PC 


+  2  A* 


(A. 12) 
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Figure  captions 


Fig.  1.  Relation  between  incident  and  scattered  phonons  in  two-dimensions 

Fig.  2.  Roton-roton  scattering  process  in  two-dimensions. 

Fig.  3.  Temperature  variation  of  the  roton  part  of  *(T). 

Fig.  4.  Temperature  variation  of  the  various  characteristic  times. 

Fig.  5.  The  coefficient  k (T)  of  thermal  conduction  as  a  function  of 
temperatures. 

Fig.  6.  The  coefficient  h (T)  of  first  viscosity  vs  temperature. 


